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ABSTRACT 



The interpretation of Fourier spectra in the time domain is critically 
examined. Power density spectra defined and calculated in the time domain 
are compared with Fourier spectra in the frequency domain for three different 
types of variability: periodic signals, Markov processes and random shots. The 
power density spectra for a sample of neutron stars and black hole binaries are 
analyzed in both the time and the frequency domains. For broadband noise, the 
two kinds of power spectrum in accreting neutron stars are usually consistent 
with each other, but the time domain power spectra for black hole candidates 
are significantly higher than corresponding Fourier spectra in the high frequency 
range (10-1000 Hz). Comparing the two kinds of power density spectra may 
help to probe the intrinsic nature of timing phenomena in compact objects. 



Subject headings: Stars: binaries — Stars: neutron — X-rays: stars — Methods: 
data analysis 



INTRODUCTION 



The Fourier transform is widely used in timing analysis in astronomy. Through a 
discrete Fourier transform, a light curve x(tk) can be decomposed into sine wave components 
in the frequency domain 

a(f j ) = J2x(t h )e~ i2 ^ kAt , (1) 
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where a(fj), j = —N/2, • • • , N/2 — 1, is the Fourier amplitude of sine wave components at 
frequency fj = jAf with Af = l/(NAt) and x(tk) represents the photon counts during a 
time interval (tk, 4 + At), k — 1, 2, • • • , N, with t^ — {k — l)At. The Fourier spectrum can 
be a powerful technique in the search for periodic signals, pulsations, and quasi-periodic 
oscillations in the light curve of a X-ray source. The Fourier power density 

Pf(fj) = Wj)\ 2 (2) 

is used to describe the variability amplitude at different frequency fj. A peak with power in 
excess of the noise distribution power in the Fourier spectrum indicates the existence of a 
periodic component in the light curve. When estimating the amplitude of a periodic signal 
from the Fourier power density spectrum, one must assume that the shape of a signal pulse 
is sinusoidal, which is usually not true for variability due to real physical processes. 

Even more care should be taken in interpreting the Fourier spectrum of an aperiodic 
process in the time domain. It is a prevalent misconception that the Fourier power spectrum 
is the only way to express the distribution of variability amplitude vs. time scale, and that 
Pf(fj) is a quantity describing the variability amplitude for the time scale 1/ fj. The actual 
relation between the Fourier power and the variability process occurring in the time domain 
is given by Parseval's theorem: 

1 N/2-1 N 

T7 E Wi)l 2 = EI*(**)l 2 > ( 3 ) 

JV j=-2/N k=l 

or equivalently, 

EP/(/i)A/ = £|r(f*)| 2 Af , (4) 

j k 

where r(t^) = x(tk)/ At is the counting rate at tk- Parseval's theorem states that the integral 
of the Fourier power density over the whole frequency range is equal to the variability 
power of the same process in the time domain. In fact, that is just the reason why the 
quantity Pf(fj) = \a>(fj)\ 2 is called power density. However, Parseval's theorem says nothing 
about the power density distribution over the time domain. The Fourier power at any 
given frequency within the limited range, — > I /At, actually represents a sum of complex 
amplitudes from an infinite number of aliased frequencies, each contribution of which is 
reduced by a factor of sinc(7r/At), with the sum of the squares of all the factors being 
unity. In contrast, the power at any given frequency in a time-based spectrum is most 
sensitive to the rate of change of the Fourier spectrum with frequency near the Nyquist 
frequency (/ = 1/(2 At)) and its aliases, since the sine function given above is changing 
most rapidly for such frequencies. The rms variation vs. time scale of a time series may, in 
fact, differ substantially from its Fourier spectrum. Even in the simplest case, where the 
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signal is purely sinusoidal with a frequency /, the Fourier spectrum is a 5 function for the 
continuous Fourier transform or approximately a 5 function for the discrete case. Such a 
function has little power density at frequencies except those near /, which, although in this 
case is a clear representation of what physical process is going on, is still not an accurate 
representation of the variability of the amplitude in a time series. Only when the time 
step is set equal to or much greater than the period, At > 1/f, can the corresponding 
light curve (or pulse profile) of a sinusoidal signal be completely flattened to make the time 
variation vanish. At time scales shorter than 1/f, however, a variation of intensity still 
exists in the light curve, and under some circumstances one would wish that power over this 
region should not be zero. In fact, different but mathematically equivalent representations 
with different bases or functional coordinates in the frequency domain do exist for certain 
light curves. The Fourier transform with its trigonometric basis is just one of many possible 
transforms between the time and frequency domains, and may not necessarily represent the 
true power distribution of a physical process in the time domain. 

One can also argue that any observable physical process always occurs in the time 
domain, and this is why real variability amplitudes at different time scales are useful in 
understanding a time varying process, in addition to the understanding gleaned from the 
conventional Fourier spectrum. An algorithm to calculate the power density spectrum 
directly in the time domain without using the Fourier transform has been proposed ( |Ci| 
[200ip . We introduce the algorithm and compare the power spectra in the time domain 
with the Fourier spectra for different kinds of time process model in §1. We have studied 
the power density spectra of a sample of neutron stars and black hole binaries, in both 
the frequency and the time domains. The results are presented in §3. A discussion of the 
potential gain in the understanding the intrinsic nature of physical processes occurring near 
compact objects via comparison of the two kinds of power density spectrum is given in §4. 



where the light curve x(k), k = 1, . . . , N, is a counting series obtained from a time 
history of observed photons with a time step At, r(k) is the corresponding rate series, and 
x = Y.k=i x(k)/N and f = J2k=i r (k)/N are the average counts and count rates, respectively. 
The unit used for a power P in spectral analysis in astronomy is usually rms 2 , where 'rms' 
refers to the analysed time series x. In the case of x being a counting series, the power P 



2. 



POWER SPECTRA IN THE TIME DOMAIN 



The initial definition of the variation power in a light curve x(k) is 
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can be expressed in units of counts 2 /s 2 . -P(At) represents an integral of variability power 
for the region of time scale > At. It is easy to prove that P(Ati) > P(At 2 ) if At 2 > Ati. 
Thus we can define the power density p(At) in the time domain as the rate of change of 
P(At) with respect to the time step At 

in counts 2 /s 3 or rms 2 /s. With Eq. (5) and (6), we can calculate the power density spectrum 
in the time domain for a light curve. In practice, the differential calculus in Eq. (6) can be 
performed numerically 

Fy ' At 2 -Ah w 
with P(Ati) and P(At 2 ) being the two powers at the time scale Ati and At 2 (At 2 > Ati), 
respectively, and At = (At x + At 2 )/2. 

To detect the source signal in a power spectrum against a noise background, we need 
to know the noise power P no i se (At) of a time series consisting only of noise. For a white 
noise series where the x(k) follows the Poisson distribution, the noise power is 

P (At) - Var(x) - <X> - 1- (R) 
nmse{ ' (At) 2 (At) 2 At' { ' 

where < x > is the expectation value of x, but with units of counts 2 , which also the variance 
of the Poisson variable x, and r is the expectation value of the counting rate which can 
be estimated by the global average of counting rate of the studied observation. The noise 
power density at At = (At x + At 2 )/2 is 

/ a ,\ Pnoise(Ati) _P rao j se (At 2 ) ^ T , 

p noise {lXt) - — _ — - V) 

The signal power density can be defined as 

Vsignal(At) = P (At) - Pnolse (At) (10) 

and the fractional signal power density as 

p'^(At) = P-^ At \ (11) 

in units of 1/s or (rms/mean) 2 /s, where both 'rms' and 'mean' refer to the time series. 

To study the signal power density in the time domain over a background of noise in 
an observed photon series, we divide the observation into M segments. For each segment 
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i, the signal power density Pi, S ignai(At) is calculated by Eq. (10). The average power 
density of the studied observation is p s = Y^=\P i, signal I M and its standard deviation is 

a (Ps) = \jYjiL\{Pi, signal — Ps) 2 I '(M(M — 1)). We can use the statistical methods based on 
the normal distribution to make statistical inference, e.g. significance test, on p s . At short 
time steps, the number x of counts per bin may be too small for it to behave as a normally 
distributed variable. But it is easy to get a large enough total number (M) of segments 
from a certain observation period to satisfy the condition for applying the central limit 
theorem as well as for using the normal statistics on the mean p s . 

To compare the power spectrum in the time domain with the Fourier spectrum for the 
same process, we study three different kinds of time series. 

(1) Periodic signal 

We use Eq.(5) and (6) to calculate the power density spectrum in the time domain for 
a sinusoidal process with a period of 0.5 s. In Figure 1, a piece of the counting rate curve 
of the studied process is shown in the top panel and the power density distribution of time 
scale of the sinusoidal signal is shown by the solid line in the middle panel. From the rate 
curve, we derive a corresponding counting series with a time step 2 ms and use 8192-point 
FFTs to get the Fourier power densities Pf(fj) and the corresponding densities in the time 
area pf(Atj = 1/ fj) = Pf(fj)ff- The Fourier spectrum is shown by the dashed line in the 
middle panel of Fig. I. The areas under the two power density curves in Fig. 1 are the 
same, as dictated by Parseval's theorem. It is obvious that the Fourier spectrum can not 
be interpreted as the distribution of variability amplitude vs. time scale. The power being 
concentrated at the sinusoid frequency / does not mean that the intensity variation of the 
process exists only at the corresponding time scale 1//. From the light curves with time 
steps At = 0.05 s and 0.28 s, shown in the bottom panel of Fig. 1, we can see that variation 
of the intensity definitely exists at time scales At < 0.5 s, but almost no power does in the 
Fourier spectrum at those time scales. In contrast with the Fourier spectrum, the power 
density spectrum determined in the time domain (the solid line in the middle panel of 
Fig.l) gives a proper description of the variability amplitude distribution of time scales. 

Figure 2 shows the power density spectrum (the solid line in the bottom panel) derived 
in the time domain and the corresponding Fourier spectrum (the dashed line in the bottom 
panel) for a periodic triangular signal (the top panel). The Fourier spectrum is obtained 
with 8192-point FFT for a light curve with time step 1 ms. The many peaks at short time 
scales in the Fourier spectrum (high frequency harmonics) do not mean that there really 
exist strong variations in the process at those time scales, which are just necessitated by 
mathematically decomposing the triangular signal into sine waves. As a result, a Fourier 
spectrum may overestimate power densities at short time scales (high frequency harmonics 



- 6- 



region) for a periodic signal with pulse shape being far from a sinusoid. 
(2) Stochastic Shots 

We now consider the power spectra of a signal s(t) consisting of stochastic shots. 
Both the rise and decay fronts of the shots are an exponential with a time constant r 
taken uniformly from the range between 5 ms and 0.2 s. The separation between two 
successive shots is exponentially distributed with an average separation d = 0.5 s. The 
average signal rate is assumed as f s = 300 cts/s. The peak height of the shot follows a 
uniform distribution between zero and the maximum. The top panel of Fig. 3 shows a 
piece of the signal counting curve with time step At = 0.01 s and the solid line in the 
bottom panel of Fig. 3 is the power density spectrum in the time domain expected for 
the signal s(t). A simulated 2000 s light curve x{k) with 1 ms time step is produced by 
a random sampling of the signal curve with Poisson fluctuations plus a white noise with 
mean rate 5000 cts/s. The middle panel of Fig. 3 shows a piece of the light curve with 
At = 0.01 s obtained from the simulated 1 ms light curve. In calculating the signal power 
density at a time scale At, The total light curve with time step At was divided into M 
data segments with iV = 100 bins each. If the segments number M < 100 in cases of large 
time scales, we let M = 100 and decreased the number N of data points in each segment 
accordingly. For each segment of the light curve with time bin At\, the total powers at two 
time scales Ati and At 2 = 2Ati through Eq. (5) were calculated. The corresponding noise 
powers were calculated by Eq. (8) with r set equal to the average counting rate. The total 
power density p(At) and noise power density p no i se (At) at At = 1.5 At i can be calculated 
by Eq. (7) and (9), respectively. Finally, the signal power density in the time domain is 
Psignai(At) = p(At) — p noise (At) . In the bottom panel of Fig. 3, the plus signs mark the 
average signal power densities at different time scales. For the same light curve with 1 ms 
time binning, we also calculated the Leahy density w(fj) = 2\aj\ 2 /a for each T = 4.096 s 
segment, where aj is the Fourier amplitude at frequency fj = j/T determined from a 
4096-point FFT. It is well known that the noise Leahy density w no i se = 2, so the signal 
Leahy density can be written as w S i gna i(fj) = w(fj) — 2 and the Fourier power density of the 
signal is expressed by PF,signai(fj) — w signa i(fj)a /T ([Leahy et al 1983] , |van der Klis 1988|) . 
The dashed line in the bottom panel of Fig. 3 shows the average Fourier power density 
spectrum of the signal with respect to time scale Atj = 1/ fj. As the characteristic time 
r of a shot is taken in the range between 5 ms and 0.2 s, there should exist considerable 
variability over this time scale range as represented by the power density spectrum in the 
time domain expected for the signal s(t) (the solid line in the bottom panel of Fig. 3) 
and by that obtained from the simulated data including noise (pluses in the bottom panel 
of Fig. 3). Figure 4 represents the power spectra for another shot model with shorter 
characteristic time constant r between 0.5 ms and 2 ms, average separation between two 
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successive shots 3 ms, average signal rate f s = 500 cts/s, and noise rate 5000 cts/s. /,From 
Fig. 3 and Fig. 4 we can see that for a random shot series the Fourier spectrum is more or 
less consistent with the power spectrum in the time domain at time scales greater than the 
characteristic time scale of the model, but significantly underestimates the power densities 
at shorter time scales. 

(3) Markov process 

The Markov process or autoregressive process of first order can describe the character 
of the variability for many physical processes. A Markov process can be expressed as the 
following stochastic time series 

u{k) = a ■ u(k — 1) + e , 

where e is a Gaussian random variable with zero mean and unit variance, and the relaxation 
time of the process is r = —At/ log \a\ with At being the time step. The observed light 
curve for the signal is 

s(k) = c ■ u{k) + r s At , 

where r s is the average rate of the signal. We make a light curve of the signal s(k) with 
At = 0.01 s, r = 0.1 s, r s = 2000 cts/s and c = 2.5. A piece of the produced signal light 
curve is shown in the top panel of Fig. 5. The final observed light curve x(k) = s(k) + n(k) 
with n(k) being a Poisson noise with mean rate 5000 cts/s, is shown in the middle panel of 
Fig. 5. In the bottom panel of Fig. 5, the solid line shows the power density distribution 
expected for the signal s(k), the plus signs indicate excess power densities in the light curve 
x{k) estimated by Eqs. (5)-(10) and the dashed line by FFT. From Fig. 5, we can see, 
similar to the case of stochastic shot process, that the proposed algorithm of evaluating 
power densities in the time domain is capable of extracting the power spectrum of the 
signal from noisy data and that the Fourier spectrum significantly underestimate the power 
densities at time scales around or shorter than the characteristic time scale of a stochastic 
process. 



3. ACCRETING NEUTRON STARS AND BLACK HOLES 

As revealed by above study based on simulations, the power density spectrum derived 
in the time domain can describe the real power distribution with respect to time scale for 
different time processes. The two kinds of power spectrum, Fourier spectrum and spectrum 
in the time domain, differ in ways depending on models of the processes. By comparing the 
two power spectra, more information about the nature of the variability of an object can be 
extracted. For this purpose we calculate both the power density spectra in the time domain 
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and Fourier spectra for the publicly available data of the Proportional Counter Array 
(PCA) aboard the Rossi X-ray Timing Explorer (RXTE) for a sample of X-ray binaries, 7 
neutron stars and 7 black hole binaries. The PCA observations of X-ray objects included 
in our study are listed in Table 1. For each analyzed observation, we use the version 4.1 
of standard RXTE ftools to extract the PCA data. At 18 time steps between 0.001 s and 
2.5 s, we make the corresponding light curves with ~ 2000 s duration in an energy band as 
noted in Table 1. Then we remove all ineffective data points caused by failure due to the 
satellite, detector or data accumulation system, and calculate the power density spectra in 
the time domain, being based on a similar procedure in our simulation. The corresponding 
Fourier spectra are constructed by using 2~ 12 s (~ 244 fis) time resolution light curves 
divided into parts containing 8192 bins. 

Figure 6 shows the results from 6 observations for 5 neutron star binaries: KS 1731-260, 
4U 1705-44, GS 1826-24, 4U 0614+091, 4U 1608-522 in the low state and 4U 1608-522 in 
the high state. The two kinds of power spectrum in the studied accreting neutron stars 
are generally consistent with each other, at least for the continuum dominated region. 
The feature is rather complicated for such sources whose Fourier spectra have significant 
quasi-periodic oscillation (QPO) structure. The left panel of Fig. 7 is from the QPO 
source, Sco X-l, whose power density spectrum in the time domain also shows a strong 
QPO structure but the peak shifts to shorter time scales due to the steep slope of the sine 
function near the Nyquist frequencies (7 - 12.5 Hz) as more and more of the 6+ Hz QPO is 
accommodated by the time sampling. The difference between the spectrum of time domain 
and the Fourier spectrum at short time scales for Cyg X-2 (see the right panel of Fig. 7) 
may also be caused by QPOs. 

For black hole candidates, we first analyze the canonical source Cygnus X-l. On May 
10, 1996 (day 131 of 1996), the All-Sky Monitor on RXTE revealed that Cyg X-l started a 
transition from the normal low (hard) state to a high (soft) state. After reaching the high 
state, it stayed there for about 2 months before going back down to the low state (|Cui et al 



1997 ). During this period, 11 pointing observations of Cyg X-l were made by RXTE. We 
use one observation of PCA/RXTE for each of the four states of Cyg X-l. The signal power 
density spectra in the time domain of Cyg X-l are shown by the plus signs in Fig. 8 and 
the corresponding Fourier spectra by the dots in Fig. 8. The power spectra of Cyg X-l in 
different states (Fig. 8) have a common trait that the Fourier spectra are significantly lower 
than the corresponding power spectra of time domain in the time scale region of At < 0.1 s. 

Besides Cyg X-l, the black hole candidates GRS1915+105, GRO J1655-40, GX 339-4, 
and XTE J1550-564 also demonstrate a significant excess in the power spectra in the time 
domain in comparison with the corresponding Fourier spectra at time scales shorter than 
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~ 0.1 s (see Figure 9). But the another analyzed black hole binary, GRS 1758-258, with 
QPO structures in its Fourier spectrum behaves differently (Fig. 10). 

4. DISCUSSION 

A Fourier power density spectrum presented in the time domain can not be interpreted 
as the real power density distribution of the physical process studied. In principal, the 
power densities in the time domain can be derived from a Fourier spectrum only when one 
knows every power density spectrum in the time domain for each sinusoidal function at all 
Fourier frequencies and adds them up with weight factors being the Fourier amplitudes. 
We propose here studying power spectra directly in the time domain. The definition of the 
power -P(At) (Eq. 5) is based just on the original meaning of rms variation and the power 
density spectrum p(At) (Eq. 6) represents the distribution of the variability amplitude 
vs. time scale. The power density spectrum p(At) in the time domain obtained from an 
observation depends only on the intrinsic nature of the signal process and the statistical 
property of the observed data, as does the Fourier power spectrum in the frequency domain. 
Our simulation studies show that the proposed algorithm, Eqs. (7), (9) and (10), is capable 
of extracting power densities of the signal from noisy data (comparing the expected signal 
spectra, the solid lines in the bottom panels of Figs. 3 and 5 and in Fig. 4, with the 
spectra calculated from the noisy data, plus signs in corresponding plots). These results 
indicate that the technique of spectral analysis in the time domain is a useful tool in timing 
and worth applying in temporal analysis for different sources. From Figs. (1) - (5) one 
can see that the difference of the Fourier spectrum with the spectrum in the time domain 
is dependent on the model of time series and sensitive to the characteristic time of a 
stochastic process. We can then use the difference between two kinds of spectrum to study 
the intrinsic nature of a studied process. 

The power spectra shown in Figs. (8) and (9) for black hole candidates have a common 
character that the Fourier spectra are significantly lower than the corresponding time-based 
spectra in the time scale region of At < 0.1 s, which is similar to the simulation results 
for the stochastic shot model (Fig. 3) or the Markov process (Fig. 5). The existence of 
stochastic shots in X-ray light curves of Cyg X-l has been noticed for a long time. With 
an improved searching and superposing algorithm and PCA/RXTE data of Cyg X-l in 
different states, Feng, Li & Chen (1998) find that the average shot profiles can be described 
by exponentials with characteristic time scales ~ 0.01 — 0.1 s. Pottschmidt et al (1998) 
point out that an autoregressive process of first order with a relaxation time r of about 0.1 s 
can reproduce approximately the variability of Cyg X-l. Thus, the characteristics of time 
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process in Cyg X-l revealed from our power spectral analysis is consistent with that from 
modeling its light curves. The black hole candidates GRS 1915+105, GRO J1655-40, GX 
339-4 and XTE J1550-564 also demonstrate characteristics similar to Cyg X-l, indicating 
that a stochastic process with a characteristic time ~ 0.1 s may be common in accreting 
black holes. While the absence of the broad-band noise above approximately 100 Hz in 
black hole candidates has been noticed before QSunyaev fc Revnivtsev 2000|) , our analysis 



shows that the apparent absence in Fourier spectra is caused by the existence of a stochastic 
process with characteristic time ~ 0.1 s and by the insensitivity of the Fourier technique to 
detecting rapid variability in a stochastic process. 

At the same time, the two kinds of power spectrum are more or less consistent for 
accreting neutron stars with continuum dominated Fourier spectrum, as shown in Fig. 6. 
Our simulation results, Figs. (3) - (5), show that for a stochastic process the two kinds of 
spectrum can be consistent with each other at time scales greater than the characteristic 
time of the process. Assuming that a significant variability of an accreting neutron star 
comes from a stochastic process with very short characteristic time constant r < 1 ms 
can explain the consistence of the two kinds of power spectrum observed in neutron star 
systems. Sunyaev & Revnivtsev (2000) find that the power density spectra of accreting 
neutron stars with a weak magnetic field have significant broad noise component at the 
frequency 500-1000 Hz. They suggest that those X-ray transients which demonstrate 
significant noise in their X-ray flux at frequencies above ~ 500 Hz should be considered 
neutron stars. Most sources studied in this work (neutron star systems 4U 0614+091, 4U 
1608-522, GS 1826-24, 4U 1705-44, KS 1731-260, Cyg X-2, and black hole candidates Cyg 
X-l, GX 339-4, GRS 1915+105, GRO J1655-40, XTE J1748-288, GRS 1758-258) are also 
studied by Sunyaev & Revnivtsev (2000) and our results support their claim under the 
condition that the spectra and noise in their statement are restricted within the Fourier 
framework. In contrast to the results from Sunyaev & Revnivtsev (2000), our results are 
based on inferring the characteristic time from the relation between two kinds of spectrum, 
no matter what the absolute magnitude of power spectrum is. The characteristic feature we 
find appears in all spectral states of the black hole candidate Cyg X-l and the same is true 
with the neutron star binary 4U 1608-522. 

The characteristic features in rapid variability for two different kinds of X-ray binary 
in our sample, which associate significant stochastic processes at the time scale ~ 0.1 s 
for black hole candidates and <C 1 ms for neutron stars, are revealed only in continuum 
or broad noise. The QPO components behave in a more complicated fashion, but can be 
understood through the sensitivity of the time-based spectrum to variation at 2St, i.e., at 
half the Fourier sampling frequency (the Nyquist frequency). It has been found that QPO 
structures in power spectra can be caused by different kinds of signal, e.g., modulated 
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periodic signals and stochastic autoregressive processes with order > 2 (|van der Kris 1986|) . 
In the case of the neutron star binary Sco X-l (the left plot of Fig. 7), the power density 
spectrum in the time domain reveals the QPO feature surprisingly well, though with a peak 
shifting to shorter time scales and with worse resolution. 

Different processes with essentially different natures could result in almost the same 
Fourier power spectrum, thus distinguishing them is difficult through timing analysis only 
with the Fourier technique. This is a reason why we need to develop and apply alternative 
methods to supplement the Fourier technique in spectral analysis. Our results, though only 
preliminary, show that simultaneous use of both the Fourier and the time domain methods 
can help in probing the intrinsic nature of timing phenomena, and further, in distinguishing 
between different kinds of accreting compact object. 

The authors thank the referee for helpful comments and suggestions and Dr. Qu Jinlu 
for help in data treatment. This work is supported by the Special Funds for Major State 
Basic Research Projects and the National Natural Science Foundation of China. The data 
analyzed in this work are obtained through the HEASARC online service provided by the 
NASA/GSFC. 
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Table 1. The used PCA observations of X-ray binaries 



Type 


Object 


Observation 


Band (keV) 


Note 




KS 1731-260 


10416-01-02-00 


3-21 






4U 1705-44 


20073-04-01-00 


3-20 






GS 1826-24 


30054-04-01-00 


3-20 




Neutron Star 


4U0614+091 


30054-01-01-01 


3-21 






4U1608-522 


30062-01-01-04 


3-21 


low state 




4U1608-522 


30062-02-01-00 


3-21 


high state 




Sco X-l 


30035-01-02-000 


2-18 






Cyg X-2 


30418-01-01-00 


2-21 






Cyg X-l 


10412-01-01-00 


2-13 


low to high 




Cyg X-l 


10512-01-08-00 


2-13 


high state 




Cyg X-l 


10412-01-05-00 


2-13 


high to low 




Cyg X-l 


10236-01-01-03 


2-13 


low state 


Black Hole 


GRS1915+105 


20402-01-05-00 


5-22 






GRO J1655-40 


20402-02-25-00 


5-22 






GX 339-4 


20181-01-01-00 


4-22 






XTE J1550-564 


30191-01-14-00 


2-13 






GRS 1758-258 


30149-01-01-00 


3-21 






XTE J1748-288 


30185-01-01-00 


2-21 
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Fig. 1. — Distribution of power density vs. time scale of a sine signal. Top panel: counting 
rate curve of the sine signal with period 0.5 s. Middle panel: power densities. Solid line - 
power density spectrum derived in the time domain, dashed line - Fourier spectrum, derived 
by 8192-point FFT for the sine signal light curve with step 2 ms and shifted downwards by 
a factor of 100. Bottom panel: three light curves with time step 0.05 s, 0.28 s, and 0.5 s, 
respectively. 
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Fig. 2. — Distribution of power density vs. time scale of a periodic triangle signal. Top 
panel: counting rate curve of the periodic triangular signal with period 0.5 s. Bottom panel: 
power densities. Solid line - power density spectrum derived in the time domain, dashed 
line - Fourier spectrum, derived by 8192-point FFT for the the light curve of the signal with 
1 ms step and shifted downwards by a factor of 80. 
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Fig. 3. — Distribution of power density vs. time scale of a shot model. Top panel: the 
signal, stochastic exponential shots with time constant r between 5 ms and 0.2 s. Middle 
panel: simulated data which involved both signal and Poisson noises. Bottom panel: signal 
power densities . Solid line - power density distribution of time scale expected for the signal. 
Dashed line - excess Fourier spectrum from the simulated data. Plus signs - excess power 
densities calculated by the timing technique in the time domain for the simulated data. 
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Fig. 4. — Distribution of power density vs. time scale of a shot model. Time constants 
of exponential shots are between 0.5 ms and 2 ms. Solid line - a theoretical distribution 
of power densities expected for the signal. Dashed line - the excess Fourier spectrum from 
the simulated data. Plus - the excess power densities calculated in the time domain for the 
simulated data. 
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Fig. 5. — Distribution of power density vs. time scale of a Markov process. Top panel: the 
signal, a Markov process with relaxation time 0.1 s. Middle panel: simulated data which 
include signal and Poisson noises. Bottom panel: power densities of the signal. Solid line 
~ the power distribution of time scale expected for the signal. Dashed line - excess Fourier 
spectrum from the simulated data. Plus - excess power densities calculated in the time 
domain for the simulated data. 
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Fig. 6. — Power density vs. time scale of 5 accreting neutron stars. Plus: spectrum in 
time domain. Dot: Fourier spectrum. 
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Fig. 7. — Power density vs. time scale of two accreting neutron stars. Plus: spectrum in the 
time domain. Dot: Fourier spectrum. 
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Fig. 8. — Power density vs. time scale of Cyg X-l. Plus: spectrum in the time domain. Dot: 
Fourier spectrum. 
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Fig. 9. — Power density vs. time scale of 4 accreting black holes. Plus: spectrum in the time 
domain. Dot: Fourier spectrum. 
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Fig. 10. — Power density vs. time scale of black hole binary GRS1758-258. Plus: spectrum 
in the time domain. Dot: Fourier spectrum. 



